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We calculate moments of the 0(al) heavy flavor contributions to the Wilson coefficients 
of the structure function F2(x,Q 2 ) in the region Q 2 S> m 2 . The massive Wilson coef- 
ficients are obtained as convolutions of massive operator matrix elements (OMEs) and 
the known light flavor Wilson coefficients. The calculation of moments of the massive 
OMEs involves a first independent recalculation of moments of the fermionic contribu- 
tions to all 3-loop anomalous dimensions of the unpolarized twist-2 local composite 
operators stemming from the light -cone expansion [JJ. 

1 Introduction 

Deep- inelastic scattering (DIS) of charged or neutral leptons off proton and deuteron targets, 
in the region of large enough values of the gauge boson virtuality Q 2 = —q 2 , allows to 
measure the leading twist parton densities of the nucleon, the QCD-scale Aqcd, and the 
strong coupling constant a s (Q 2 ) = a s (Q 2 )/(47r), to high precision. For unpolarized DIS 
via single photon exchange, the double-differential cross section can be expressed in terms 
of two inclusive structure functions p2,i(x, Q ). These decompose for twist r = 2 into a 
Mcllin convolution of non-perturbative masslcss parton densities fj(x, /J 2 ) and the perturbative 
Wilson coefficients Cj( 2 ,L)( x , Q 2 /^ 2 , m k/^ 2 )- The latter describe the hard scattering of the 
photon with a massless parton. They are given by the sum of the purely light - denoted by 
Cj,(2,L) ~~ an d heavy flavor contributions, Hj ( 2 .l). Here k — c, b and j — q, g, depending on 
the type of process one considers, x denotes the Bjorken scaling variable. Especially in the 
region of smaller values of Bjorken-cc, the structure functions contain large cc-contributions 
of up to 20-40 %, denoted by F^ c h {x^ Q 2 )- The perturbative heavy flavor Wilson coefficients 
corresponding to these structure functions are known at NLO semi-analytically in x-space 
[2]. Due to the size of the heavy flavor corrections, it is necessary to extend the description 
of these contributions to O(o? s ), and thus to the same level which has been reached for the 
massless Wilson coefficients [3]. 

A calculation of these quantities in the whole kinematic range at NNLO seems to be out 
of reach at present. However, in the limit of large virtualities Q 2 , Q 2 > 10 m 2 in the case 
of F£ c (x,Q 2 ), one observes that F£ c L (x,Q 2 ) are very well described by their asymptotic 
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expressions [1] neglecting power corrections in m 2 /Q 2 . In this kinematic range, one can 
calculate the heavy flavor Wilson coefficients analytically. This has been done for FSf 'ix, Q 2 ) 
to 2-loop order in [4j |5] and for Ff^(x : Q 2 ) to 3-loop order in [6]. Note that in the latter 
case, the asymptotic result becomes valid only at much higher values of Q 2 . The asymptotic 
expressions are obtained by a factorization of the heavy quark Wilson coefficients into a 
Mcllin convolution of massive OMEs Ajk and the massless Wilson coefficients Cj t i, if one 
heavy quark flavor of mass m and rif light flavors are considered. In the present paper, we 
report on the calculation of the massive OMEs Ajk to 3-loop order for fixed even moments 
of the Mellin variable N, cf. [7] for details. We further calculate the OMEs which are 
required to define heavy quark parton densities in the variable flavor number scheme [S]. 
We also obtain moments of the terms oc Tp of the 3-loop unpolarized anomalous dimensions 
jij . Our results agree with those obtained in [3] . Since the present calculation is completely 
independent by method, formalism, and codes, it provides a strong check on the previous 
results. 

2 Heavy Flavor Operator Matrix Elements 

The heavy flavor Wilson coefficients for a single massive quark may be expressed as 
uS.ps.ns ( Q 2 n% 2 \ _ ^g.ps ( 0^ , r s,PS,NS ( 0^ rrf\ 

where the photon couples to a light (L) or heavy (H) quark line, respectively. Further S 
stands for the flavor-singlet contributions, which are separated into a pure-singlet (PS) and 
non-singlet (NS) part via S = PS + NS. The factorization formula for the inclusive Wilson 
coefficients reads in Mellin space, [H [H] > 

^S.PS.NS.as/,, Q 2 m 2 \ s-^ ,S,PS,NS/, r ™ 2 \ ^S.PS.NS / AT . -. Q 2 \ /„n 

Here /z refers to the factorization scale between the heavy and light contributions in Cjj2,L) 
and 'as' denotes the limit Q 2 m 2 . The Cji2,L) are precisely the light Wilson coefficients 
and describe all the process dependence. The arguments (n/), (n/ + 1), indicate at how 
many light flavors the respective quantities have to be taken. This factorization is only valid 
if the heavy quark coefficient functions are defined in such a way that all radiative correc- 
tions containing heavy quark loops are included. Otherwise ([2]) would not show the correct 
asymptotic Q 2 -behavior [S]. The mass dependence is given by the process independent 
massive OMEs Aij , which are the flavor-decomposed twist-2 operator matrix elements 

2 °° 2 

' 1=1 ^ 

Here, i denotes the external on-shell particle [i = q, g) and Ok stands for the quarkonic 
(fc = q) or gluonic (fc = g) operator emerging in the light-cone expansion. The subscript H 
indicates that we require the presence of heavy quarks of one type with mass m. The logarith- 
mic terms in m 2 / fi 2 are completely determined by renormalization and contain contributions 
of the anomalous dimensions of the twist-2 operators. Thus at NNLO the fermionic parts 
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of the 3-loop anomalous dimensions calculated in Refs. [5] appear. All pole terms of the 
unrenormalized results provide a check on our calculation and the single pole terms allow 
for a first independent calculation of the terms oc Tp of the 3-loop anomalous dimensions. 

In case of the gluon operator, the contributing terms are denoted by A gq> Q and A gg ^Q. 
For the quark operator, one distinguishes whether the operator couples to a heavy or light 
quark. In the NS-case, the operator, by definition, couples to the light quark. Thus there is 
only one term, A^q. In the S and PS-casc, two OMEs can be distinguished, {A^q, A S q} 
and {^4qq, A^ g }, where, in the former case, the operator couples to a light quark and in 
the latter case to a heavy quark. 

Eq. ^ allows to calculate the heavy flavor Wilson coefficients in the limit Q 2 ^> m 2 up 
to 0(al) by combining the results obtained in Ref. [5] for the light flavor Wilson coefficients 
with the 3-loop massive OMEs which are computed in this work [Jj. 

A related application of the heavy OMEs is given when using a variable flavor number 
scheme to describe parton densities including massive quarks. The OMEs are then the 
transition functions going from rif to rif + 1 flavors. One thus may define parton densities 
for massive quarks, see e.g. Ref. [S]. This is of particular interest for heavy quark induced 
processes at the LHC, such as cs — > W + at large enough scales Q 2 . 



3 Renormalization 



We work in Feynman gauge and use dimensional regularization in D = 4 + e dimensions, 
applying the MS-scheme, if not stated otherwise. The renormalization proceeds in four steps, 
which we will briefly sketch here and refer to [7] for more details. Mass renormalization is 
performed in the on-shell scheme [10] , whereas for charge renormalization we use the MS- 
scheme. We work in an intermediate MOM -scheme for charge renormalization by requiring 
that the heavy quark loop contributions to the gluon propagator vanish for on-shell external 
momentum. This is necessary for the renormalization of the massive OMEs to cancel infrared 
singularities which would otherwise remain. Z g in this MOM-scheme can be calculated using 
the background field method [IT] . Finally, we transform our result back to the MS-scheme 
for coupling constant renormalization via, [7J, 
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The remaining singularities are of the ultraviolet and collinear type. The former are renor- 
malized via the operator Z-factors, whereas the latter are removed via mass factorization 
through the transition functions T. After coupling- and mass renormalization, the renor- 
malized heavy flavor OMEs are then obtained by 



A = Z- 1 AT~ 1 



(6) 



where quantities with a hat are unrenormalized. Note that in the singlet case operator 
mixing occurs and hence Eq. ^ should be read as a matrix equation, contrary to the NS- 
case. The Z— and T-factors can be expressed in terms of the anomalous dimensions of the 
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twist-2 operators to all orders in the strong coupling constant, cf. pJ[T2] up to O(a^). From 
Eq. ([5]) one can infer that for operator renormalization and mass factorization at O(a^), 
the anomalous dimensions up to NNLO [9] together with the 1-loop massive OMEs up to 
0(e 2 ) and the 2-loop massive OMEs up to 0(e) are needed. The 2-loop OMEs up to O(e ) 
were calculated in Refs. [4j (5] [13]. Higher orders in e enter since they multiply Z— and 
T-factors containing poles in e. This has been worked out in detail in Ref. [12], where we 
presented the 0(e) terms a,Q g , a^g 5 and ag,^- The terms a g 2 g q and a^Q were given in 
Refs. |13j . Thus all terms needed for the renormalization at 3-loops in the unpolarizcd case 
are known. 

Finally we would like to point out the difference between the MOM- and MS-scheme for 
coupling constant renormalization at NLO [13]. Eq. ([2]) holds only for completely inclusive 
quantities, including radiative corrections containing heavy quark loops [8]. Additionally, 
([2]) has to be applied in such a way that renormalization of the coupling constant is carried 
out in the same scheme for all quantities contributing, i.e., the MS-scheme. If one evaluates 
the heavy-quark Wilson coefficients, diagrams of the type shown in Fig. [T] may appear 
as well. It contains a virtual heavy quark loop correction to the gluon propagator in the 
initial state and contributes to the terms L g ^ and respectively, depending on whether 
a light or heavy quark pair is produced in the final state. Note that in the former case, 
this diagram contributes to F^ 2 .l) { x i Q 2 ) ' ln t ne inclusive case, but is absent in the semi- 
inclusive QQ-production cross section. In Refs. [2], the coupling constant was renormalizcd 
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Figure 1: 0(a 2 s ) virtual heavy quark corrections to hr 2 
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in the MOM-schcmc in 0(a 2 s ) by absorbing the contributions of the above diagram into the 
coupling constant. This can be made explicit by considering the complete gluonic Wilson 
coefficient up to 0{a 2 s ), including one heavy quark, see Eq. @, 



C g ,2{n f ) + L g - 2 {n f + 1) + H g , 2 (n f + 1) = a 
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The above equation is given in the MS-scheme. Here, the diagram shown in Fig. Q] con- 
tributes, corresponding exactly to the color factor Tp. Transformation to the MOM-scheme 
for a s , Eq. JH, yields 



C fll2 (n/) + L gt2 (n f + 1) + H gt2 (n f + 1) 



4 (l),MOM n^lr, _L 1 "I 
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In the above equation, all contributions due to diagram [T] have canceled, i.e. the color factor 
Tp does not occur at the 2-loop level in the MOM-scheme. Splitting up Eq. (jHJ) into H g ^ 
and Lg t i, one observes that L g ^ vanishes at 0(a 2 s ). The term H g i is the one calculated in 
Ref. [4], which is the asymptotic expression of the gluonic heavy flavor Wilson coefficient 
as calculated exactly in Refs. [2]. It is not clear whether the same can be achieved at the 
3-loop level as well, i.e., transforming the general inclusive factorization formula ((2]) in such 
a way that only the contributions due to heavy flavors in the final state remain. Therefore 
one should use the asymptotic expressions at 3 loops only for completely inclusive analyzes. 
This approach has also been adopted in Ref. [H] for the renormalization of the massive 
OMEs, which was performed in the MS-scheme and not in the MOM-schcmc, as previously 
in Ref. [4]. In the NS-case a similar argument holds, which can be found in Ref. [4]. 



4 Calculation and Results 

The massive OMEs at O(a^) are given by 3-loop self-energy type diagrams, which contain a 
local operator insertion. The external massless particles are on-shell. The heavy quark mass 
sets the scale and the spin of the local operator is given by the Mellin-variable N. The steps 
for the calculation are the following: We use QGRAF 14J for the generation of diagrams. 
Approximately 2700 diagrams contribute to all the OMEs. For the calculation of the color 
factors we refer to |15j . The diagrams are then genuinely given as tensor integrals. Applying 
a suitable projector provides the results for the specific Mcllin moment under consideration. 
The diagrams are further translated into a form, which is suitable for the program MATAD 
[16j . through which the expansion in e is performed and the corresponding massive three- 
loop tadpole-type diagrams are calculated. We have implemented all these steps into a 
FORM -program [17] and checked our procedures against various complete two-loop results 
and certain scalar 3-loop integrals and found full agreement. 

Applying Eq. (O, one can predict the pole structure of the unrenormalized results and 
thus the logarithmic terms of the renormalizcd OMEs. These contributions can be expressed 
in terms of the anomalous dimensions up to 3 loops, the expansion coefficients of the QCD 
/3-function up to 2 loops and the 1- and 2-loop contributions to the massive OMEs. Thus 
the logarithmic terms are known for general values of N. This is not the case for the constant 

(3) 

term, which contains the genuine 3-loop contributions a\j . These are known for the fixed 
values of N as calculated in this work [7J. 

For the OMEs A ( Q ] g , A q 3 g ] Q and A { gg \ Q the moments N = 2 to 10, for ^q]' PS to N = 12, 

(3) NS ('3') PS (3) 

and for A qq Q , A qq q , A gq Q to N = 14 were computed. For the flavor non-singlet terms, 
we calculated as well the odd moments N = 1 to 13, corresponding to the light flavor 
— combinations. The complete calculation took about 250 days of computer time. All our 
results agree with the predictions obtained from renormalization, providing us with a strong 
check on our calculation. As an example, we show the constant term of A gq ^Q for N = 2 
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Here Q denotes the Riemann ^-function at integer argument i and the term is given by 



It appears in all OMEs we calculated and is known to arise as a genuine mass effect. 
5 Conclusions and Outlook 

We calculated all massive 3-loop OMEs for even Mellin-moments TV = 2. ..10(12, 14) using 
MATAD. This confirms for the first time, in an independent calculation, the moments of 
the fermionic parts of the corresponding 3-loop anomalous dimensions [5]. Combining our 
results with [3] , this provides fixed moments of the heavy flavor Wilson coefficients of F2 in 
the limit Q 2 3> to 2 . First phenomcnological studies of the effects of our calculation are in 
preparation. 
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